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Abstract: In this paper we describe a formalism for generating exclusive final states in 
diffractive excitation, based on the optical analogy where diffraction is fully determined by 
the absorption into inelastic channels. The formalism is based on the Good- Walker formal- 
ism for diffractive excitation, and it is assumed that the virtual parton cascades represent 
the diffractive eigenstates defined by a definite absorption amplitude. We emphasize that, 
although diffractive excitation is basically a quantum-mechanical phenomenon with strong 
interference effects, it is possible to calculate the different interfering components to the 
amplitude in an event generator, add them and thus calculate the reaction cross section 
for exclusive diffractive final states. The formalism is implemented in the DIPSY event 
generator, introducing no tunable parameters beyond what has been determined previously 
in studies of non-diffractive events. Some early results for DIS and proton-proton collisions 
are presented, and compared to experimental data. 
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1. Introduction 

Although most events in high energy colhsions have a continuous distribution of hadrons, 
events with rapidity gaps are not rare; they contribute of the order 10% both in DIS at 
Hera [1,2], and in pp coUisions at the CERN SppS cohider [3,4] and the Tevatron [5]. 
These events have commonly been interpreted as the shadow of absorption into inelastic 
channels, in analogy with diffraction in optics. In this view diffraction is expected to be 
a quantum-mechanical phenomenon, to be analyzed at amplitude level, and not directly 
treatable by semiclassical probabilistic methods. 

In the Regge formalism the absorption is represented by cut pomeron or reggeon di- 
agrams. Regge diagrams, cut in between two exchanged pomerons or reggeons, represent 
elastic scattering, while diffractive excitation originates from triple pomeron diagrams, cut 
through one of the pomerons [6,7]. At high energies multi-Regge diagrams, where the rela- 
tive contributions from cut and uncut diagrams are given by the AGK cutting rules [8] , are 
important. Recent analyses within this formalism include work by Kaidalov and cowork- 
ers [9], Ostapchenko [10], the Durham [11,12], and the Tel Aviv groups [13-15]. These 
analyses give results for inclusive diffractive cross sections and distributions in da/dM-^, 
but do not give information about details in the diffractive final states. (Ref. [9] makes 
the extra assumption that a pomeron interacts like a qq pair.) A scheme for generation of 
exclusive diffractive final states is outlined in ref. [16], but as far as we know, no results are 
presented yet. We note that, in contrast to our approach which is based on amplitudes, in 
this scheme diffractive events are generated as a semi-classical probabilistic process. 

Diffractive excitation has also been described within the Good- Walker formalism [17], 
where it is the result of differences in absorption probability between different components 
of the projectile wavefunction. In the Regge approaches mentioned above, this formalism 
was used only for low mass excitation. It was, however, early proposed by Miettinnen and 
Pumplin [18], that fluctuations between parton cascades with different absorption probabil- 
ity also can describe excitation to high masses. In QCD a high energy proton is visualized 
as a virtual cascade, with partons filling the rapidity range between the proton rest frame 
and the observer. This implies that the proton wavefunction also contains components 
with high masses ~ exp{rapidity range). In high energy collisions the cascades are de- 
termined by the BFKL dynamics, which has a stochastic nature with large fluctuations, 
and in the dipole formulation these cascades are also eigenstates of the interaction [19]. 
These fluctuations were studied by Hatta et al. [20] in an analysis of diffractive excitation 
in DIS at very high energies, within the saturated region where < Q% This approach to 
diffractive excitation has also been exploited by the Lund group within the dipole cascade 
model implemented in the DIPSY Monte Carlo [21,22], with applications to DIS and pp 
collisions at present day collider energies. It is the aim of this paper to generalize this 
approach to calculate amplitudes for transitions to exclusive final states. 
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In both the triple-Regge and the Good- Walker formahsms the difFractive ampUtude is 
via the optical theorem determined by the inelastic cross section. At high energies inelastic 
events result from gluon exchange, which causes colour connections between projectile and 
target, and the diffractive amplitude is represented by the (uncut) perturbative BFKL 
pomeron formed by a two-gluon ladder. The stochastic nature of the BFKL pomeron is 
also present in the AGK cutting rules. In ref. [23] it is argued that the triple-pomeron and 
Good- Walker formalisms only represent different views of the same phenomenon. This idea 
is also supported by the fact that the bare pomeron in the DIPSY model also reproduces 
the triple-Regge form for the diffractive excitation cross section [22]. 

In the "colour reconnection" approach it is assumed that the colour exchange from an 
initial hard subcollision can be neutralized through the mediation of subsequent soft glu- 
ons [24,25]. In this picture the cross section for gap events is determined by a reconnection 
probability. This can be tuned to fit experimental data, but is in this picture not dynam- 
ically fixed by the total inelastic cross section, and the relation to diffraction in optics is 
not clear. 

Of particular interest in analyses of diffractive excitation are the increased effects 
of saturation at higher energies. In DIS diffractive excitation of the photon has been 
calculated from clastic scattering of qq and qqg states. Golec-Bicrnat and Wiisthoff [26,27] 
have emphasized that the coupling of a virtual photon is significant to large qq dipolcs with 
highly asymmetric energies, but that the coupling to the proton is suppressed for dipolcs 
larger than the saturation scale i?o ~ ^/Qs{x). This implies on one hand a suppression of 
the diffractive cross section, and on the other hand that the interaction is not dominated 
by soft interaction, but can be treated by perturbative QCD. This argument should also 
be applicable for diffraction in pp collisions, suggesting that also this can be treated by 
perturbative methods. 

In pp collisions the absorptive cross section is larger than in DIS. In the Good- Walker 
formalism diffractive excitation is determined by the fluctuations in absorption probability, 
which become small when the black disk limit is approached. Thus diffraction is dominated 
by elastic scattering, when the saturation effects become large at higher energy, in partic- 
ular for small impact parameters. In the Regge formalism these saturation effects are 
represented by "enhanced diagrams" , which interfere destructively and reduce the proba- 
bility for a rapidity gap. Also in the colour reconnection formalism, the large number of 
partons participating in the collision tends to fill in a potential gap. In the analysis by 
Goulianos [28], these saturation effects are described by a "renormalized pomeron flux", 
which starts to become noticeable in pp collisions around -^i k. 20 GeV. 

Most of the analyses cited above estimate inclusive distributions in da/dM'j^ (or 
da / drjrasLx) or semi-inclusive hard processes. To understand the properties of exclusive final 
states additional information or assumptions are needed. At lower masses, Mx, diffr ac- 
tively excited pions or protons fragment in a string-like manner; an excited pion is similar 
to an e"^e~" -annihilation event [29], while an excited proton has similarities to DIS [30,31]. 
These observations are in agreement with the model by Donnachie and Landshoff, in which 
the pomeron is assumed to interact with the quarks in the target like a photon [32]. 

For higher masses, the limited acceptance of the detectors at the CERN SppS col- 
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lider and the Tcvatron implies that the experimental data usually do not cover the full 
phase space, and the experimental knowledge about the properties of high mass diffractive 
final states is therefore also limited. At the CERN SppS collider rapidity distributions 
were measured by the UA5 collaboration [4] . The result was consistent with a p^-limited 
fragmentation, although there was here no experimental information about transverse mo- 
menta. Similar results were obtained by the UA4 collaboration [3], which pointed out the 
similarity between a diffractive system and a non-diffractive pp collision at ^/s = Mx- At 
the Hera detectors the coverage of the photon fragmentation end is quite good, although 
an excited proton is mainly outside the acceptance. Here the diffr actively excited photon 
gives higher multiplicity and lower thrust than an e'''e~-annihilation event with the same 
mass [33]. 

High p± jets in diffractive states were observed by the UA8 collaboration [34], but the 
limited acceptance did not allow a detailed study of the diffractive state, e.g. of how the 
recoil was distributed. "Hard diffraction" with high p± jets have later also been observed 
at the Tevatron [35,36] and at Hera [37,38]. Due to the implications for detection of 
exclusive central Higgs production, the observation of hard central production with two 
rapidity gaps at the Tevatron, e.g. of dijets [39] and e'^e~ pairs [40], has gained special 
interest. 

Hard diffraction has frequently been analyzed within the Ingclman-Schlein formalism 
[41], assuming that the pomcron interacts as composed of partons, in a way similar to a 
hadron. The diffractive interaction is expressed in terms of a flux factor fwpixp), describing 
the probability to find a pomeron with energy fraction xp in the proton, times the inelastic 
interaction between the pomcron and the projectile proton or virtual photon. In this 
way analyses of Hera data for both hard and soft diffraction have been described by 
a common set of "pomeron structure functions" /g/p,p(/3), which have been fitted to a 
DGLAP evolution scheme [42]. These analyses indicate a dominant gluonic component in 
the pomeron. This leads to more complicated final states than the straight strings seen at 
lower masses, in agreement with the observations in refs. [3] and [33] mentioned above. 

This approach is also implemented in several Monte Carlo (MC) event generators used 
in analyses of experimental data, e.g. Pompyt [43], Rapgap [44], and PythiaS [45]. In 
these models the result depends on pomeron flux factors and pomeron structure functions, 
which have to be fitted to data. Here the effects of saturation have to be included in the 
flux factor, or in "gap survival factors" . 

The aim of this paper is to present a formalism, in which the cross section for exclusive 
diffractive final states can be calculated directly based on the dynamical features of small x 
evolution and saturation. The formalism is an extension of the Lund Dipole Cascade model 
implemented in the DIPSY MC [46,47], which is based on BFKL evolution, including es- 
sential non-leading effects, and saturation within the cascade evolution. Hence the model 
is based on perturbative QCD, which might be motivated by Golec-Biernat-Wiisthoff's 
argument that gluons below the saturation scale Qg are suppressed at small x. Diffractive 
excitation is calculated within the Good- Walker formalism with no additional free param- 
eters. Results for inclusive diffractive cross sections have been presented previously [21,22], 
the generation of exclusive non-diffractive final states was discussed in ref. [47], and the 
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model is here extended to describe exclusive final states in diffractive excitation. 

We want to emphasize that, although in the optical analogy diffractive excitation is 
basically a quantum-mechanical phenomenon with strong interference effects, it is in our 
formalism possible to calculate the different interfering contributions to the amplitude in 
the DIPSY MC, add them with their proper signs, and thus calculate the reaction cross 
section. This distinguishes our formulation from other approaches, in which final states 
are generated as a semi-classical probabilistic process. 

The outline of the paper is as follows: In section 2 we review shortly the eikonal 
approximation and the Good-Walker formalism. The dipole cascade model is described in 
section 3, and in section 4 we discuss how this model can be extended to describe exclusive 
final states in diffractive excitation. In section 5 we discuss the implementation in the 
DIPSY event generator, and in section 6 we present some early results. Some possible 
future developments are discussed in section 7, followed by our conclusions in section 8. 

2. Diffractive excitation in the Good— Walker picture 
2.1 The Good— Walker formalism 

In the Good-Walker formalism [17], diffraction in hadronic collisions is analogous to diffrac- 
tion in optics, and a result of absorption and unitarity. Unitarity constraints and saturation 
are most easily treated by the eikonal formalism in impact parameter space. 

Assume that the elastic scattering is driven by absorption into a number of inelastic 
states n, with Born amplitudes \/2fn- In the eikonal approximation [48] the absorption 
probability, Pabs) is given by the inelastic cross section 

Pabs = dair^Jcfb = 1 - where F = Y,fn- (2.1) 

Thus the probability for the projectile not to be absorbed into any of the inelastic states 
is given by e~^^. 

For a structureless projectile unitarity then gives the S'-matrix 

S = e-^ = e-^f". (2.2) 

Thus the amplitude for elastic scattering is given by 

T = 1 - 5 = 1 - e~-^, (2.3) 

which satisfies the optical theorem T = ^(T^ + Pabs)- Note that we have here defined the 
amplitude T without the conventional factor i, so that T becomes real. The elastic cross 
section is given by dae\/d^h = T^, and adding the inelastic cross section, the total cross 
section is given by datot/d^b = 2(1 — e~^) = 2T. 

If the projectile has an internal structure, it can be excited in a diffractive scattering 
event. This implies that the mass eigenstates (asymptotic incoming and outgoing states) 
can differ from the diffractive eigenstates, i.e. eigenstates to the 5-matrix. We denote the 
diffractive eigenstates with eigenvalues given by S^i = e~^»$i. They can be absorbed 
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with probability 1 — e ^ % and have elastic scattering amplitudes (<I>j|T|<I>j) = 6ijTi, with 

The diffractive eigenstates can be written as linear combinations of the mass eigen- 
states, 

<^i = J2^ik'i'k, (2.4) 

k 

where Cik is a unitary matrix (thus = Y1 '^L*^«)' ^^'^ incoming state is given by = 
^Q. The elastic amplitude is then given by the average over the diffractive eigenstates: 

(^'o|r|M/o) = 5^|c,o|'ri = (r), (2.5) 

with the elastic cross section 

daei/(fb = Iciol'Ti)' = {Tf. (2.6) 
The amplitude for diffractive transition to the mass eigenstate becomes 

(^-fclTl^o) = J]4TiQo, (2.7) 

i 

which gives a total diffractive cross section (including elastic scattering) 

dadm/dH = Y,{^o\T\^k){'fk\T\^o) = (T^)- (2.8) 

k 

Subtracting the elastic scattering we find the cross section for diffractive excitation 

da^m eJdH = da^m/dH - da,,/dH = {T^) - {Tf = Vt, (2.9) 
which thus is determined by the fluctuations in the scattering process. 
2.2 Diffractive eigenstates 

The basic assumption in our earlier analyses of diffractive excitation in refs. [21,22], is that 
the diffractive eigenstates correspond to parton cascades, which can come on shell through 
interaction with the target. As mentioned in the introduction, this was also the assumption 
in the early work by Miettinen and Pumplin [18], and a similar approach has been used by 
Hatta et al. [20]. This means that the factors qo in eq. ( p.4| ) are given by the evolution of 
the virtual cascade. The process is illustrated in fig. |l|. Fig. ||a shows the virtual cascade 
before the collision, and fig. illustrates an inelastic interaction, where gluon exchange 
gives a colour connection between the projectile and the target. This implies that the 
beam is absorbed with probability 1 — e~'^^\ In this and the following figures, solid lines 
represent real emissions present in the final state, while a dashed line corresponds to a 
virtual emission, which did not come on shell via the interaction with the target, and 
therefore is reabsorbed in the cascade. Fig. ||c shows an elastic interaction. Via the optical 
theorem the elastic amplitude is given by the inelastic cross section, represented by the 
square of the diagram in fig. ||b, which corresponds to the exchange of two gluons. Note 
that the elastic amplitude in eq. ( |2.5[ ) is the result of coherent interaction of all partons in 
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Virtual cascade 
a 



Inelastic int. 
b 



Elastic scatt. 
c 



Diffractive ex. 
d 



Figure 1: (a) An example of a parton (or dipole) cascade evolved in rapidity, (b) The exchange 
of a gluon gives rise to an inelastic interaction, (c) Elastic scattering is obtained from coherent 
scattering of different partons in different cascades, via the exchange of two gluons. (d) Diffractive 
excitation is obtained when the result of the two-gluon exchange does not correspond to the coherent 
initial proton state. Dashed lines indicate virtual emissions, which are not present in the final state. 



all possible cascades. Fig. |T](i, finally, shows the contribution to the scattered beam, which 
is orthogonal to the incoming state. This corresponds to diffractive excitation, with the 
amplitude given by eq. (2/7). The lines can symbolise gluons in a traditional cascade, or 
dipoles in a dipole cascade. When the diagram in fig. is squared, it is consistent with 
Mueller's triple Regge formalism [6]. 



3. The dipole cascade model 
3.1 Mueller's dipole model 

Mueller's dipole cascade model [49-51] is a formulation of LL BFKL evolution in transverse 
coordinate space. Gluon radiation from the colour charge in a parent quark or gluon is 
screened by the accompanying anticharge in the colour dipole. This suppresses emissions 
at large transverse separation, which corresponds to the suppression of small k_i_ in BFKL. 
For a dipole with charges in transverse points x and y, the probability per unit rapidity 
(Y) for emission of a gluon at transverse position z is given by 

dV a 2 (x - y)^ ., , _ NcUs , . 

— = —d z- — -, with Q = . 3.1) 

dV 2-K (x — zj^(z — yj^ TT 

As shown in fig. |2| the dipole is split into two dipoles, which (in the large Nc limit) emit 
new gluons independently. The result is a cascade, giving a dipole chain where the number 
of dipoles grows exponentially with Y. 

When two cascades collide, a pair of dipoles with coordinates (xj,yj) and (xj,yj), in 
the projectile and target respectively, can interact via gluon exchange with the probability 
2fij, where 



fij = /(xi,yi|xj,yj) 



/(x,-y,)2(y^2xn2 



log 



(Xj -Xj)2(y, -yj^2 



(3.2) 



We note here that the interaction probability goes to zero for a small dipole. This implies 
that the singularity in the production probability for small dipoles in eq. (|3.1| ) does not 
give infinite cross sections. We note also that gluon exchange means exchange of colour 
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Figure 2: Gluon emission splits the dipole into two dipoles. Repeated emissions give a cascade, 
which produces a chain of dipoles. 




Figure 3: An interaction between a dipole in the projectile and another in the target due to gluon 
exchange gives a recoupling of the dipole chains. 




y 



Figure 4: Double interaction results in a dipole loop, corresponding to a pomeron loop. 

between the two cascades. This implies a reconnection of the dipole chains, as shown in 
fig. ^, and the formation of dipole chains connecting the projectile and target remnants. 

In Mueller's model the constraints from unitarity are satisfied using the eikonal for- 
malism. When more than one pair of dipoles interact, colour loops are formed, as shown 
in fig. I 

3.2 The Lund dipole cascade model DIPSY 

It is well known that non-leading effects are very important in BFKL evolution. Another 
problem is that Mueller's model does not include pomeron loops within the cascade evo- 
lution. The Lund model [46, 47, 52, 53] is a generalization of Mueller's model, which also 
includes NLL BFKL effects, nonlinear (colour-suppressed) effects in the evolution, and 
confinement effects. 

a) Beyond LL BFKL 

The NLL corrections to BFKL evolution have three major sources [54]: 



-8- 




Figure 5: Two dipoles with tlie same colour form a colour octet, which may be better approximated 
by dipoles formed by the closet colour-anticolour pairs. This implies a recoupling of the dipole 
chains, which favours the formation of small dipoles. Thus the number of dipoles resolved by a 
probe with given resolution is reduced. 

Non-singular terms in the splitting function: These terms suppress large z-values in 
the individual parton branchings. Most of this effect is taken care of by including energy- 
momentum conservation. This is effectively taken into account by associating a dipole with 
transverse size r with a transverse momentum k± = 1/r, and demanding conservation of 
the light-cone momentum p-|- in every step in the evolution. This gives an effective cutoff 
for small dipoles. 

Projectile-target symmetry: A parton chain should look the same if generated from the 
target end as from the projectile end. The corresponding corrections are also called energy 
scale terms, and are essentially equivalent to the so called consistency constraint [55]. 
This effect is taken into account by conservation of the negative light-cone momentum 
components, p-. 

The running coupling: Following ref. [56], the scale in the running coupling is taken 
as the largest transverse momentum in the vertex, 
b) Nonlinear effects in the evolution 

As mentioned above, multiple interactions produce loops of dipole chains correspond- 
ing to pomeron loops. Mueller's model includes all loops cut in the particular Lorentz 
frame used for the analysis, but not loops contained within the evolution of the individual 
projectile and target cascades. As for dipole scattering the probability for such loops is 
given by as, and therefore formally colour suppressed compared to dipole splitting, which 
is proportional to q = NcOLs/tt- These loops are therefore related to the probability that 
two dipoles have the same colour. Two dipoles with the same colour form a quadrupole 
field. Such a field may be better approximated by two dipoles formed by the closest 
colour-anticolour charges. This corresponds to a recoupling of the colour dipole chains. 
The process is illustrated in Fig. ^, and we call it a dipole "swing". With a weight for 
the swing which favours small dipoles, we obtain an almost frame independent result. The 
number of dipoles in the cascade is not reduced, and the saturation effect is a consequence 
of the smaller interaction probability for the smaller dipoles. Thus the number of dipoles 
(or gluons) resolved by a probe with a given resolution Q^, is reduced. In this way the 
swing also generates effectively 2 —t- 1, or in some cases a 2 — t- 0, transitions. 
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c) Confinement 

As mentioned earlier, saturation effects imply that the cascade evolution is dominated 
by relatively small dipoles. However, although suppressed, the rare large dipoles generated 
in a purely perturbative evolution with massless gluons give non-negligible effects, and 
eventually Proissart's bound will be violated [57]. Therefore confinement is also important, 
and is taken into account by giving the gluon an effective mass. 

3.3 Inclusive cross sections and non-difFractive final states 

The Lund cascade model is implemented in a Monte Carlo event generator called DIPSY, 
with applications to collisions between electrons, protons, and nuclei. An incoming virtual 
photon is here treated as a qq pair, with an initial state wavefunction determined by QED. 
For an incoming proton we make an ansatz in form of an equilateral triangle of dipoles, 
but after evolution the result is rather insensitive to the exact form of the initial state. 

The model reproduces well the inclusive total and elastic cross sections in pp collisions 
and DIS, as described in ref. [53]. Diffractive cross sections are reproduced within the 
Good- Walker formalism [22]. It should, however, here be noted that when the nonlinear 
effects from saturation are switched off, the result agrees with the triple-pomeron formula 
for a bare pomeron with intercept a(0) = 1.21 and an approximately constant triple- 
pomeron coupling g^p = 0.3 GeV~^. This result is a consequence of the common underlying 
assumption, with diffraction as the shadow of absorption. 

The model is also generalized to describe the properties of exclusive non-diffractive final 
states [47]. We here note that BFKL evolution properly reproduces inclusive observables. 
For exclusive final states it is necessary to take into account colour coherence and angular 
ordering, as well as soft radiation, related to the z = 1 singularity in the gluon splitting 
function. These effects are taken into account in the CCFM formalism [58,59], which also 
reproduces the BFKL result for inclusive cross sections. 

An essential point is here the fact that the softer emissions in the CCFM formalism can 
be resummed, and the total cross section, as well as the structure of the final states, are fully 
determined by the " fcx-changing" gluons. We denote the real emitted gluons in a ladder 
qi and the virtual links fcj, and momentum conservation then implies kxj_i = k^j -|- q^j. 
For fcj^-changing emissions k±i is either much larger or much smaller than k±i^i. This 
also means that q±i ~ max(k±i,k±i^i). These emissions are called "primary gluons" in 
ref. [60], and "backbone gluons" in ref. [61], and the weight for such a backbone chain is 
given by 

weight = Yl a-^dvi. (3.3) 

As discussed in ref. [62], this feature also gives a dynamical cutoff for small q±, which grows 
slowly with energy, and gives a dynamical description of the cutoff for hard subcoUisions 
needed in event generators like Pythia [45] or Herwig [63]. 

To generate the inelastic final states, we thus have to go through the following steps: 

• Generate two dipole cascades for the projectile and the target. 
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• Determine which dipoles in the projectile and target become colour connected via 
gluon exchange. 

• Extract the backbone chains from the cascades, and remove virtual branches, which 
cannot come on shell by the interaction. 

• Add softer emissions as final state radiation, with appropriate Sudakov form factors. 

• Hadronize the gluonic chains. 

The results from MC simulations, presented in ref. [47], are able to give a fair reproduction 
of data for minimum bias and underlying events from LHC and the Tevatron, remembering 
that this implementation does not include matrix element corrections for hard scattering 
or contributions from quarks. 



4. Exclusive final states in diffractive excitation 
4.1 Basic formalism 

We will now discuss how to describe exclusive states in diffractive excitation, within the 



dipole cascade model and the Good- Walker formalism. As discussed in sec. 2.2, the par 



tonic cascades are interpreted as the diffractive eigenstates. Thus the unitary matrix Cj^ in 



eq. (2.4) is represented by a unitary operator, describing the cascade evolution. The cas- 
cade is absorbed with probability 1 — e~^^, and via the optical theorem this gives an elastic 
scattering amplitude T = 1 — . The hermitian conjugate matrix, c|,^, in the amplitude 
for diffractive transition to another mass eigenstate in eq. (2.7), represents an evolution 



backwards in rapidity, which can reabsorb some of the emitted partons. As cascades with 
more partons generally are absorbed with higher probability, this backward evolution does 
not always give the original proton back, but a different partonic system, interpreted as a 
possible mass eigenstate. 

The evolution in rapidity may also be interpreted as an evolution in time, from an 
incoming state at t = — oo to the cascade present just before the collision at time t = 0. 
Thus the states at t = — oo correspond to the mass eigenstates in sec. (we include 
the possibility that this state does not correspond to the ground state ^'o)) while the 
evolved states at time t = —6 correspond to the diffractive eigenstates ^i. As the possible 
cascades are not a discrete set, we describe the cascade evolution by a unitary operator 
U{—6, — oo), replacing the unitary matrix Cik in eq. ( |2.4f ). The interaction of the cascades 
with the target is described by a diagonal operator T, which corresponds to an evolution 
operator U{+5, —6) = f/jnt = I — T (where I denotes the unit matrix). The transformation 
back to mass eigenstates, obtained by the inverse of the cascade evolution operator, can 
then be interpreted as describing the evolution from time t = +6 to t = oo, U{+oo, +6) = 
U{—6, — oo)"'^ = U{—6, — oo)''". Thus the total S'-matrix is given by 

S = U{+oo, +6)U{+6, -6)U{-6, -oo) = U{-6, -oo)-'^UintU{-6, -oo), (4.1) 
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with the diffractive final state given by 

^-out = S%n. (4.2) 

In order to simphfy the subsequent equations, we here expressed the relations in terms of 
the S'-matrix, instead of the transition amplitude T. Note that as the diffractive scattering 



is the shadow of absorption, this operator in not unitary. In section 7A an extended 
formalism is discussed, in which the non-diffractive states are included, and the 5'-matrix 
is a unitary operator. 

In order to illustrate the essential features in diffractive excitation, we first study some 
simple toy models. The results will help understanding the more realistic gluon cascades 
discussed in the later subsections. In the toy models in sec. |4^, and in the continuous 



cascade in sec. 4.3, the gluons are assumed to interact individually, with no effects of 



screening from neighbouring gluons. Such screening effects are taken into account in the 



dipole cascade discussed in sec. We will here also assume that the cascade only includes 
gluon emission, and that effects of saturation from gluons joining (sometimes referred to 
as gluon recombination) can be neglected. Such effects are discussed in sees. ^ and 7^. 



4.2 Toy models 

4.2.1 A two-particle system 

We first study the simple case of a particle, the valence particle, which can emit a single 
gluon at a fixed point in phase space. (This example is also directly analogous to the 
2-channel Good-Walker treatment of low mass excitation in the work by the Tel Aviv 
or Durham groups.) The state with no emission and only the original valence particle 
is denoted |1,0), and the state including the emitted gluon is denoted Thus the 

numbers indicate whether the two particles are present (1) or absent (0). At asymptotic 
times t = iboo these states are also mass eigenstates, |l,0)±oo and |l,l)±co respectively. 
During the evolution up to the collision, the particle can emit a gluon. This is described 
by the unitary evolution operator U{—5,—oo), which contains a |1,0) — )■ |1,1) transition. 
In the (|1, 0), |1, 1)) space it can be described by a unitary matrix (we choose phases such 
that U is real): 

U{-6,-^)=(j^~^y a2 + /32 = i. (4.3) 

Thus gives the probability for gluon emission, and the probability for no emission. 

In the collision, which takes place in a short time interval {—5, +S), we assume that the 
valence particle and the gluon can be absorbed with probabilities 1 — e"^-'^^ and 1 — e~^-^^ 
respectively. In this interaction the beam particle becomes colour-conected to the target, 
and thus absorbed from the incoming beam. This implies that the states |1,0) and |1, 1) 
have probabilities e~'^-^° and e~'^^-^''^-^'^^ not to be absorbed. In accordance with the opti- 
cal theorem, the evolution operator describing the evolution of the diffractive eigenstates 



during the collision, is given by (c/ eq. (|2.2|) ): 



U^n, = U{+5, -5) (4.4) 
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Figure 6: The diagrams for elastic scattering in a two particle system. In this and subsequent 
figures a dashed line represents a virtual emission, which is reabsorbed by its parent. A dot rep- 
resents a possible emission or absorption, which did not occur. The long dashed line indicates the 
interaction with the target, with the value for [/jnt indicated. 

Finally the system evolves from time +(5 to -|-(X), described by the hermitian conjugate 
of the unitary matrix U{—5, —00): 



U{+oo, +5) = U^{-5, -00) = " ^ (4.5) 




The S-matrix for the diffractive states is thus given by 

S = U{+oo, +5)UintU{-5, -00) 



e-^o(a2 + /32e-^i), -a/3 e-^o(l - e^^^) 
-Q/3e-^o(l - e"-''!), e-^o(/32 + Q2e"-/'i) 



(4.6) 



The scattering matrix is given by T = I — S*, and from the above expression we can read off 
the scattering amplitudes for an incoming single valence particle. Thus the amplitudes for 
elastic scattering and diffractive excitation can be written as (recalling that + = 1) 

Elastic scattering : r(|l,0) ^ |1,0)) = 1 - Su = a'^{l - e'^") + /3^{l - e-^°-^^)i4.7) 
Diffractive excit. : r(|l, 0) ^ |1, 1)) = -S'21 = a/3 e-^''(l - e'^i). (4.8) 



We note here in particular, that these results agree with the expressions in eqs. ( p.6| , 
|2.9|) , for the case with Cik = U{—6, —00) given by eq. ( ^^ ) and Tj = 1 — (C/int)^ (no 
summation) given by eq. ([4.4| ). 

It will be helpful to represent the different contributions to the amplitudes by associated 
diagrams. Thus the elastic amplitude is interpreted as the sum of the two diagrams in fig. |6|. 
In the first diagram no gluon is emitted before the collision (weight a), followed by the 
weight for no absorption of the valence particle (weight e~-^°), and finally no emission is 
allowed after the interaction (weight a). In the second diagram a virtual gluon is emitted 
(weight /?), followed by the weight for no absorption of the two-particle system (weight 
5-/0-/1)^ and finally the virtual gluon has to be reabsorbed by the valence particle (weight 

In the same way the amplitude for diffractive excitation is represented by the diagrams 
in fig. 1^. Here the emitted gluon appears as a real particle in the final state. In the first 
diagram the gluon is emitted before interaction with weight /3, then the system avoids 
absorption by the target with weight e~-^°~-^^, and finally the gluon also avoids reabsorption 
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Figure 7: The diagrams for diffractive excitation in a two particle system. 



by its parent with weight a. In the second diagram the gluon avoids emission before 
interaction with weight a, the system is not absorbed by the target with weight e"-^" , and 
finahy the gluon is emitted after the interaction with weight (— This possibility is a 
consequence of the unitarity of the evolution operator U{—6,—oo) in eq. (|4.3|) , and the 
minus sign comes because it is now emitted by and not by U. The two contributions 
interfere destructively, which gives the factor (1 — e"-^^), which also can be interpreted as 
the amplitude for elastic scattering of the emitted gluon. 

This simple example illustrates two essential features of the amplitudes for exclusive 
states in diffractive excitation: 

1. Fluctuations in the scattering process are necessary for diffractive excitation, as 



seen in eq. (2^). In this toy model the probability for gluon emission is equal to /3 , so 
the average number of emitted gluons hitting the target is (n) = • 1, and the variance is 
(n^) - (n)2 = /32 . i2 _ (^2 . ^ Q,2^2^ rpj^^g ^Yie factor a/3 in the amplitude in eq. (U) 
(to be squared in the cross section) describes the fluctuations in the emission process. The 
last factor is the difference between the interaction amplitudes of the states with and with- 



out emission, and eq. (4^) shows that diffractive excitation requires fluctuations between 
different cascades with different interaction amplitudes. We also see that diffractive exci- 
tation vanishes in the black disk limit, where /o and fi become very large, and diffractive 
scattering becomes purely elastic. 

2. Diffraction is fundamentally a quantum mechanical process, where interference is 
important. The two contributions in fig. ^ or |^ have to be added in the amplitude, and 
not in the cross section. 

It should also be noted that the time t can be replaced by any variable, which 
parametrizes the evolution between the incoming mass eigenstates and the eigenstates 
for the interaction. A BFKL based cascade would use the rapidity y, while a DGLAP 
cascade would prefer the transverse momentum k±. We will here keep the notation t, and 
change to y in sec. where we will apply the model in DIPSY. 

4.2.2 Cascade with two possible emissions 

We now study a system in which two different gluons can be emitted. The gluons are 
denoted 1 and 2, and there are then four different states, with one, two or three particles: 

|100), |110), |101), |111). (4.9) 

The two numbers are here occupation numbers for the three possible particle states. 
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We first study a sequential cascade, where the first gluon can be emitted from the 
valence particle, and the second gluon can be emitted from the first one. Thus the emissions 
are assumed to be ordered in the evolution variable, and the second gluon can only be 
emitted if the first one already is present. Note that we here do not consider the possibility 
of a gluon is emitted from one parton and then reabsorbed by another. Such recombination 
effects beyond the large Nc limit are discussed in sees. H and 7.2. 



The evolution operator has now two components: the first, Ui, allows the first gluon 
to be emitted or absorbed by the valence particle, and the second, U2, allows the second 
gluon to be emitted or reabsorbed by gluon 1. In the basis in eq. ( |4.9D , these operators can 
be written 



U = U2U1, where Ui 



( ai 

/3i Qi 

ai 

V /3i 



\ 


-/3i 
ai / 



U2 



/l \ 

02 -/32 
10 
\ /32 a2 / 



Thus the first gluon can be emitted with probability /3^, and the second with probability 
/3|, once the first is present. This asymmetry also implies that the two matrices U\ and {72 
do not commute. 

The three particles have absorption probabilities 1 — e"^-^" , 1 — e^^-^^ , and 1 — e^^-^^ 
respectively. The evolution operator describing the interaction with the target is then given 
by 



/ e-/o 



V 













-/()— /2 











-h 



\ 



(4.10) 



As in the earlier example, the evolution after the interaction is given by the inverse of 
U2U1. Thus the projection on diffractive states of the S-matrix is 



S 



ululUiMUi. 



(4.11) 



Multiplying the matrices given above, and looking at the first column in the matrix T 
I — S, we get the amplitudes relevant for an incoming single valence particle: 



T(|100) 
T(|100) 
T(|100) 



Tel = T(|100) ^ |100)) = a'iil - e-^°) + l3'tai{l - e^-^o^^^) 



5-/0-/1 



1110)) 
1101)) 



/0-/l-/2^ 

-^^)+/3|( 
e--^o/3i^e~^iQ2/32(l -e--^2) 



e 

e"^«ai/3i[ai(l 

-/o fl2„-/i 



-/1-/2 



;"-''oai/3ie~^ia2/32(l - e"^""). 



(4.12) 
(4.13) 
(4.14) 
(4.15) 



As in the previous example, these results can also be obtained from a sum of diagrams, 
as illustrated in figs. |8|~[Tl|. In these diagrams the emission before, or the absorption after, 
the interaction gives a factor /3, while an emission after interaction gives — /3. The absence 
of a possible emission or reabsorption gives a factor a. The interaction with the target 
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Figure 8: The diagrams for elastic scattering (|100) |100)) in an ordered three gluon system. 
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Figure 9: The diagrams for diffractive excitation of the first emission (|100) —> |110)) in an 
ordered three gluon system. 
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Figure 10: The diagrams for diffractive excitation of the second emission (|100) |101)) in an 
ordered three gluon system. 

gives a factor , where F = fi, with the sum running over the gluons present at the 
interaction. The sum of the contributions from the relevant diagrams give the S-matrix 



elements corresponding to the amplitudes in eqs. (4.12 - 4.15). 

In case of two gluons, which can be independently emitted from the valence quark, the 
excited state |110), with a single emission, can be created through the four diagrams in 
fig. which add up to 

T(|100) ^ |110)) = e--^»ai/3i(l - e-f'){al + /^fe'^^) (4 jg) 

The transition to state |111), with both gluons emitted, would in this case also have four 

(4.17) 



contributions, as shown in fig. |13|. The result would then be 

r(|000) ^ |111)) = e"^«ai/3ia2/32(l - e~^i)(l - e-^''). 
4.2.3 Interpretation 

The results in eqs. (0, PI) and ( plHilTl) can be represented by a few rules, which in 



sec. 4.3 and Appendix A are shown to hold also for a general cascade. 
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Figure 11: The diagrams for diffractive excitation of both emissions (|100) in an ordered 

three gluon system. 
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Figure 12: The diagrams for diffractive excitation of one out of two (|100) —?■ |110)) independent 
emissions. 
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Figure 13: The diagrams for diffractive excitation of two (|100) — >■ |111)) independent emissions. 



1. A final state is specified by the real partons present in the final state. It will also be 
important to separate those real emissions, which are the last gluon in their chain, 
from those, called parent gluons, which have emitted at least one other real gluon. 

The different diagrams, contributing to the amplitude, may also contain "virtual" 
gluons, which are emitted but reabsorbed in the cascade. Let Cy. denote the proba- 
bility for emission of a cascade Vi from the real parton i. The virtual cascades can 
contain more than one emission, and including no emission as representing an empty 
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cascade, conservation of probability implies that Cy-. = 1 for each i. The sum is 
over all possible virtual cascades emitted from the parton i. 

2. The contribution from a last real parton i is given by 

a,ft^4^(l-e-^»-^^0- (4.18) 
Vi 

3. The contribution from a parent real parton i is given by 

«.ft5]4^e-^--^^. (4.19) 
Vi 

Comments: 

1) All real emissions come with a factor Either they are emitted before the 
interaction (weight /3j) and not reabsorbed afterwards (weight ai), or not emitted before 
(weight Oj) and then emitted afterwards (weight — /3j). 

2) If no virtual emissions were possible, every last gluon gives a factor {l — e~-^^). This 
is the result from interference between contributions from emission before and after the 
interaction, and can also be interpreted as the amplitude for elastic scattering from the 
target. Examples are gluon 1 in eq. (^^) and fig. 0, and gluon 2 in eq. (l4l^) and fig. |ll]. 

3) When virtual emissions are possible, the factor (1 — e"-^*) for a last gluon is modified. 



An example is the amplitude T(|100) — t- |110)) in eq. ( [1.13 ), where gluon 2 can be emitted 
and reabsorbed by gluon 1, with probability Gluon 2 can only be emitted if gluon 
1 is emitted before the interaction, and therefore only affects the last term in the factor 
(1 — e~'^'), changing it to (1 — e~f^~f^). Here the suppression factor e~.^^ represents the 
weight for no absorption of gluon 2 from the target. (The cross section contains the factor 
e~'^-^^, which is the probability for no absorption.) The probability for not emitting gluon 2 
(or for emitting an empty cascade) is 02 = 1 — /^l) and in this case there is no modification 
of the weight (1 — e~'^^). Adding the two contributions can be expressed as the sum in 
eq. ( pSl) . 

4) The parent gluons cannot be emitted after interaction (they would then not be able 
to emit their daughters), and therefore come with a factor e~f% which can be interpreted 
as the weight for not being absorbed by the target. The effect of a virtual cascade V, is to 
reduce the weight for no absorption by a factor e~^^ , as expressed in eq. ( [4.191) . 

5) These results work also for the original valence particle, if it is treated as a "last" 
particle for the elastic amplitude, and as a "parent" particle in diffractive excitation, where 
it has emitted at least one real gluon. (There is no factor ao/3o, as the valence particle is 
always present.) The difference is only that the two terms in the parenthesis in eq. ( 4.1^ ) 



do not correspond to emissions before and after interaction. It is here a contribution from 
the definition T = I — S, which gives an extra term 1 in the elastic amplitude, reflecting 
the interference between the incoming and the scattered waves. With this interpretation, 
an example of the sum over virtual emissions from a parent particle is given by gluon 2 in 
eq. ( p^ and fig. |l2. 
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4.3 Continuous cascade 



We now consider a cascade were every gluon can emit further gluons at every point in the 
evolution parameter t. This corresponds to an infinite product of evolution operators U{t) 
with infinitesimal /3j(t). In this subsection the gluons are assumed to interact individually, 
with no effects of screening from neighbouring gluons. Such screening effects are taken 



into account in the dipole cascade discussed in sec. iA. We will here also assume that the 
cascade only includes gluon emission, and that effects of saturation from gluons joining 
within the cascades can be neglected. 

The infinitesimal values for /3j implies that the corresponding values for ai are very close 
to 1, but as there is an infinite set of ai, the small differences from 1 cannot be neglected. 
However, as the constraint a? -\- = 1 represents a conservation of probability, and 
therefore the a:s can be thought of as Sudakov factors, these corrections are automatically 
taken into account in a MC event generator like DIPSY. 

This probabilistic interpretation is also applicable for the emission of virtual cascades 
Vi, emitted from a real gluon i. Let v denote the individual gluons in the virtual cascade V, 
and V the gluons which could have been emitted, either from gluon i or from the gluons in 
Vi, but were not. In a general situation the probability for the gluon i to emit the cascade 
Vi is then given by 

Here is the weight for emission and reabsorption of the gluon v in Vi, while a? is the 
weight for no emission of gluon v in Vi (where V is the set of all possible emissions not 
included in Vi). As in the toy models discussed above, the sum of emission probabilities 
for all possible cascades from the parent i, adds up to one, 

E^v-Efn/5m«i)=i' (4-21) 

Vi Vi \v&Vi vGVi J 

as a result of the unitarity relationship between a and /3. In the sum over all cascades V 
we here also include the "empty cascade" Vq, with no emissions, for which Cy^ = Yl^^y. a^. 

In Appendix A we demonstrate that the rules presented in the previous section also 
are relevant for a general cascade, and the amplitude for the transition to a real state \R) 
is given by 



= n f^rar HE 4,e"^-^-^ n E (1 - e-f'-'^^ ) . (4.22) 
VreiJ / \peP Vp J \ieL y, / 

Here R is the set of all "real gluons". is a virtual cascade emitted from a gluon p in 
the set of "parent gluons", denoted P. / is a gluon in the set L of "last gluons", and Vi 
is a virtual cascade emitted from I. The sums over Vp and Vi run over all such cascades, 
including the empty cascades. Finally and Fy^ denote the sum of the relevant terms 
fi- 
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The probalastic interpretation of the factors Cy in eqs. ( [4.20| , 4.21| ) imply that this result 



can be generated by a MC event generator. The result can be interpreted as the probability 
to create the real state, times the probability that none of the parents is absorbed in the 
interaction with the target, times the probability that the last gluon of each branch in the 
cascade should interact elastically. 

In this section the gluons are assumed to interact individually, with no effects of screen- 
ing from neighbouring gluons. Such screening effects are taken into account in the dipole 
cascade discussed in the following section. We have also assumed that the target does not 



fluctuate like a parton cascade. Two colliding cascades will be discussed in sec. 4.5 



4.4 Prom gluons to dipoles 

4.4.1 Differences from the individual gluon cascade 

In the gluon cascade, discussed in the previous section, the interaction of a gluon with the 
target is not changed by the emission of daughter gluons. The gluons interact individually, 
with no effects of screening from neighbouring gluons. Such screening effects are taken 
into account in the dipole cascade, where the screening is determined by the compensating 
charge in the other end of a dipole. Thus a gluon in a small dipole has a smaller cross 
section, leading to colour transparency. When a dipole is emitting a gluon, the old dipole 
disappears, and is replaced by two new dipoles, as discussed in section ^. This is also the 
case when a virtual cascade is emitted, which implies that some of the real dipoles are no 
longer present in the interaction with the target. Therefore the interaction amplitude F 
cannot as easily be separated in one part from the real dipoles and one from the virtual 
emissions. Some details in the calculations are rather technical and therefore left for ap- 



pendix A. Here we present first a toy model in sec. 4.4.2, which illustrates the problem. 



and then a sketch of the full result in sec. 4.4.3. 



4.4.2 Toy model with single emission 

We consider a simple toy model with two states, an initial single dipole and a state where 
this is split in two daughter dipoles. Denoting the states by the occupation numbers of 
the three dipoles, the single dipole state is called |100) and the split state with two new 
dipoles is called |011). The weight for no absorption is e~-^°, e~-^^, and e~-^^ for the three 



dipoles. The diagrams contributing to diffr active excitation are illustrated in fig. 14. To 
get a |011) final state, the emission can happen before or after interaction, which gives two 
contributions to the amplitude. Thus we get 

r(|100) ^ |011)) = -5(|100) ^ |011)) = -[a/3e~^i-^2 +a(-/3)e"-^o] = a/3(e-^" - e'^^-^^) 

(4.23) 

We note here that the weight e~'^° cannot be factored out in the same way as in eq. ([4.^). 
It will turn out that the last dipole to split, that is the "second last dipole", will play a 
similar role also in the general case, discussed in the next section. 
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Figure 14: The amplitude for tlie transition |100) |011) contains two contributions. The dipole 
can spUt by gluon emission before (left) or after (right) the interaction with the target. 

4.4.3 Full dipole cascade 

First the notation, illustrated in fig. |l^, has to be slightly revised, as the interesting sets of 
dipoles are different in some aspects. As in the previous section, the set of real dipoles in 
the final state is denoted R. Each branch of a dipole cascade also has a last gluon, which 
can be emitted before or after the interaction, while the rest of the emissions all have to 



occur before the interaction. In fig. 15 the last gluons in the cascade are circled. Each last 
gluon is connected to two dipoles, formed when the last gluon was emitted. The set of 
all such dipole pairs is denoted L. (Note also that two different last pairs can never have 
a common dipole.) As in the cascade of independent gluons, the remaining real dipoles 
will be called parent dipoles (included in the set P), although this notation is not fully 
adequate. They are not parents in the sense that they have split into daughter dipoles. 
They do, however, connect two gluons, which both have been involved in further emissions. 



In the cascade with independent gluons in sec. 4.3, the last emission is the addition 
of one more gluon, and the previous gluons remain unchanged. As discussed above, in a 
dipole cascade the last emission will add two new dipoles, and remove the emitting dipole. 
The dipoles which are split in the last emission play a special role, and are called hidden 
dipoles. In diagrams where the last split occurs before the interaction, the hidden dipoles 
do not contribute to the weight for no absorption from the target, but they do contribute if 
the last emission happens afterwards. This feature is illustrated in the toy model discussed 
above. 

Virtual cascades in the dipole formalism will also remove the dipole it is emitted from, 
while adding more dipoles. Therefore the emission of virtual cascades from the "parent 
dipoles" do not give the factorizing non-interaction weight exp(— Fp — Fvp) obtained in 
eq. (I^H). Instead we get when the cascade Vp is non-empty, the weight exp(— Fv^), while 



for the empty cascade we get exp(— Fp). We therefore here change the notation, and include 
the empty cascade in the set of all virtual cascades, and define Fy-g = Fp. This definition 
implies that the non-interaction weight for the dipole p, including possible virtual emissions, 
is always given by exp(— Fy^,). 

We use this notation also when a "last dipole pair" / emits a virtual cascade Vi. Thus, 
when one or both of the dipoles in the pair / is "emitting" an empty cascade, the non- 
interaction weight exp(— Fvj) gets a contribution from the still intact dipole(s) in the pair /. 

Finally, in case the last dipole pair is produced after the interaction, virtual cascades 
can also be emitted by the hidden dipoles h. As these cascades must be reabsorbed before 
the last dipole can be emitted, they are restricted to the rapidity range y > yi, where yi is 
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Figure 15: A sample dipole cascade with rapidity on the horizontal axis and transverse space on 
the vertical axis. Full lines are the dipoles left in the final state, while dotted or dashed lines split 
into new dipoles before the final state. The circled gluons are the ones without children, and are 
the only gluons which can be emitted both before and after interaction. The emission of a childless 
gluon creates a pair of dipoles denoted k. The childless gluons are emitted by the "hidden" dipoles 
hi. The dipoles pi are the final state dipoles connected only to gluons with children. Dipoles that 
always split before interaction are marked with thin dotted lines, and the long dashed line is the 
original dipole. 



the rapidity of the last gluon. Also here the non-interaction weight exp(— i<V;^) is defined to 
equal exp(— F/j), when the cascade Vh is the empty cascade. With this notation the result 
in eq. (|4.22D , is for a dipole cascade replaced by the following expression: 



Tn={Uf^rar [UJ: c\e-^^^ 11 E E 4 (^""^^ " ^""^^ ) • (4-24) 
\reR / \peP Vp J \leL Vi Vn{y>yi) J 

Note here that there is no product over hidden dipoles, as the hidden dipole h is specified 
by the last pair I. The last factor, (e~^^'> - e"^^0> 

is a generalization of the result in the 
simple toy model. The details of the calculations are presented in appendix A. 

Note that, like for the independent gluon cascade in eq. ( 4.22| ), in the derivation of 



these results it is assumed that the virtual cascade from one emission is independent of the 
virtual cascade from another emission. Thus the result corresponds to a cascade without 
saturation effects from gluons which can join, or dipoles which swing. These effects are 
further discussed in sees. |5| and 



4.5 Target cascade 

Up to now the target has been described as a potential without internal structure. It can 
however also be described by a superposition of virtual cascades, Vt, in the same way as the 
projectile. For single diffractive excitation, where the target is scattered elastically, the sum 
over virtual cascades from the target particle should be summed over at amplitude level 
(c/ eq. ( |2.6D ). The amplitude can thus be written as an average over target configurations 
Vt, with weights Cy , where all interaction amplitudes F now depend on the target cascade. 
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Figure 16: A diagram contributing to single difFractive excitation in a collision between two 
cascades, where the target is scattered elastically. Labeling the partons in the upper projectile 
state 0, 1, and 2, and labeling the partons in the target cascade a, b, and c, with a the incoming 
parton, the amplitude corresponding to this diagram is /3i/32/36/3c exp(— i^)Q!ia2/3b/3c, with F = 

foa + fla + /2a + fob + fib + f2b + foc + flc + /2c- 

For a dipole cascade we then get the amplitude 



Vt \reR / UgP Vp I \leL Vi Vhiy>'i 



Vh{y>yi) 

(4.25) 

F'x,T is here the interaction amplitude between the virtual projectile cascade X and the 
target virtual cascade T, given as a sum Fx,t = X^i j fi,j '^^^^ dipoles z in X and j in T. 

5. Implementation in DIPSY 



The amplitude in eqs. can be used to generate diffractive events with the DIPSY 



event generator [47]. The approach will be to first select a real cascade with the weight 
rire-R /^^Q^r) it will appear in the squared amplitude, and then perform the sum over 
virtual cascades to determine the amplitude for this real state. 

Real Cascade. First note that in a continuous cascade, each /3 is small, and each indi- 
vidual a is very close to 1. So the weight for the real cascades is essentially Ylrenl^r- This 
weight does not include the probability that no further cascade is emitted after R. To get 
the correct weight, any time the Monte Carlo generates a certain cascade R it should count 
towards the probability for that final state, no matter if the cascade goes on to generate 
more gluons afterwards. This is done by first generating a cascade D with DIPSY as normal 
up to a maximum rapidity Umax, and return a random sub-cascade R C D, weighted by 
the total number of sub-cascades of D. 

As discussed in sec. |3.3| , the structure of inelastic final states is determined by the 
"backbone" or "/c^-changing" emissions. Presuming that diffraction is the shadow of ab- 
sorption to inelastic states, we here conjecture that also difFractive states are specified 
by their A:^-changing backbone gluons. Therefore real sub-cascades which contain softer. 
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not /c_L-changing, emissions have to go through the same reweighting procedure as the 
non-diffractive final states, as discussed in sec. 3.3 and fully described in ref. [47]. 



The sum over virtual cascades. For single diffractive excitation the sums over virtual 
cascades from the final excited state R all come with weights (?, and can be calculated 
directly from DIPSY. The sum over target cascades is independent of the excited state, 
and these cascades can be pre-generated and reused for each real state to save cpu-time. 

For each real state, the subsets P and L are identified, and virtual cascades Ip, V/, and 
Vh are generated. Each of the pre-generated target cascades Vt is paired up with a virtual 
cascade, and the factors e"'^^?"^' and {e~^^h^^t — e~^'^i'^t) in the amplitude in eq. ( 4.25| ) 



are calculated and multiplied together. The average over all the virtual cascades is then 
squared, and added to the weight for the real cascade (Hr/^r) discussed above, to give the 
total weight for the real state. 

Saturation effects. In the Lund cascade model saturation effects are included within 
the cascade evolution via the dipole "swing" (see sec. in addition to the effect from 
multiple subcollisions. The result in eq. ( 4.251) is valid for a linear cascade. In the MC 



implementation, swings are included in each virtual cascade from each subset P, h, I, but 
swings between these virtual cascades are not considered in the current implementation. 
Swings are also included in the virtual cascade Vt from the target. 

In DIS, this omission is expected to be of negligible importance. First each real chain 
end I gives a small factor (1 — e'^^'-^*)^, and therefore events with two or more chain ends 
(corresponding to multiple pomeron exchange) are suppressed, unless is very small. 
Further, a large contribution to the cross section comes from real states with or 1 emission 
from the qq pair, and in those cases there are no parent emissions P, and there can be no 
swing between different subsets of the real cascade. For states with two or more gluons, 
there are corrections for the swing, but the suppression makes it small as long as the 

emissions are not too many. To further minimize the correction, the maximum rapidity for 
the generated projectile cascades, ymaxj is chosen as close to the real state as possible, to 
give the virtual cascades from the real state little room to give saturation effects. This is 
compensated by a longer evolution of the elastically scattered target proton, which includes 
the swing. 

In the case of a diffractively excited proton, saturation is a larger effect as the gluon 
density is higher and because the average dipole size is larger. Its effect is, however, 
reduced because diffractive excitation is dominated by peripheral collisions with relatively 
few interacting dipoles. For central collisions with large interaction probability, diffractive 
scattering is dominantly purely elastic. A method to calculate the amplitude taking full 
account of saturation is introduced in section but this is not implemented in the present 
MC. 

FSR and hadronisation As in the application of the DIPSY MC for generation of 
non-diffractive final states [47], final state radiation is included in the same way as for the 
Linked Dipole Chain model [64,65] in ARIADNE [66], and finally hadronisation using the 
Lund string model [67,68] implemented in PythiaS [45,69]. 
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With these steps it is possible to generate full exclusive single diffractive final states 
from the DIPSY event generator. We want to stress that no new parameters are introduced, 
which can be tuned to experimental data. All parameters were previously determined by 
total and elastic cross section as functions of ^/s, and to some extent by exclusive non- 
diffractive data and elastic reactions. Final state radiation in ARIADNE and hadronization 
in Pythia have been tuned to exclusive LEP data. 



6. Early results 

As emphasized above, the cross sections for diffractive excitation are in our formalism, 
via the optical theorem, fully determined by the projectile and target cascades, and their 
absorption into non-diffractive inelastic reactions. For DIPSY, the inelastic absorption 
is tuned to inclusive observables and pp minimum bias multiplicity in [47], and thus no 
additional parameters are introduced in the extension to diffractive final states. It should 
also be noted that tuning to experiments was done with respect to pp data only, the only 
exception being the total 7*p cross section (or equivalently F2) as function of Q^. 

In this section we show some early results with low statistics as proof of concept. We 
will continue with improvements in the present version of the MC implementation, which 
will facilitate predictions for higher collision energies and excitation masses. 

6.1 DIS 

Experimental results for diffractive excitation of the photon in DIS are available from 
HERA. Unfortunately a diffr actively excited proton is mainly outside the acceptance of 
the HERA detectors. Results for inclusive diffractive cross sections and for dasn/dM'^ 
could be easily obtained from the original DIPSY MC, and were presented in ref. [21], in 
good agreement with data from HERA. In the new formalism it is also possible to obtain 
information about the partonic content of the diffractive states. In fig. |l^ we show our 
result for the distribution in InMj. at = 120 GeV and = 24GeV^. The figure 
also shows the separate contributions from states with 0, 1 and > 2 gluons besides the qq 
pair initially coupled to the virtual photon. We see that there is a bump at Mj^ ~ Q^, 
which corresponds to final states formed by the initial qq state. This bump is followed 
by a more fiat distribution, obtained from states including also one or more gluons. As 
expected, states with more gluons become increasingly important for higher masses Mx- 
The distribution in the figure has a smooth cut off around 50 GeV, due to the Lorentz 
frame used in the simulation, which limits the rapidity range for the gluons in the excited 
state. 

Some results for the properties of the hadronic final state, obtained after final state 



radiation and hadronization, are presented in figs. ^ and 19, and compared with HERA 
data from HI [70,71]. In all cases the results are obtained for W = 120 GeV and = 
24GeV2. Infig. [l| we show the average total charged multiplicity and its fiuctuations, as 
function of excited photon mass Mx in single diffractive excitation. We see here a very 
good agreement between our results and the experimental data, both for the average and 



for the fluctuations. Fig. 19 shows the average charged multiplicity as function of rapidity. 
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Figure 17: The distribution in InM^ for single diffractive excitation in DIS a.t W = 120 GeV 
and — 24GeV^. Besides the total result, also the contributions from qq states with no giuon 
emission, 1 gluon emission, and 2 or more gluon emissions are indicated. The distribution has a 
smooth cut off above 50 GeV due to the Lorentz frame used is the simulation. 
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HERA 




Figure 18: The average total charged multiplicity (left) and its fluctuations (right) as function 
of excited photon mass Mx in single diffractive DIS at HERA [70] . 

and average energy flow as function of pseudorapidity, in two bins of Mx ■ Also here we find 
a rather good agreement with data. The fact that it agrees for both the multiplicity and 
the energy flow, indicates that also the average transverse momentum has to be correctly 
reproduced. 



6.2 pp and pp collisions 

Also in colliders for pp or pp collisions most detectors have limited acceptance for a diffrac- 
tively excited proton. At the CERN pp collider, results for distributions in pseudorapidity 
were presented from the UA5 and UA4 detectors [3,4]. Fig. shows ?7-distributions for 
single diffraction in two mass bins with {Mx) = 100 GeV and {Mx) = 140, in pp collisions 
at = 546 GeV. We see that a rapidity plateau with dN^h/drj ~ 2.5 is developed for 
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Figure 19: The average charged multiphcity as function of rapidity in bins of Mx (left) [70], and 
average energy flow as function of pseudorapidity in bins of Mx (right) [71] . The results correspond 
toW ^ 120 GeV and = 2AGcV'^ . 



(a) dNj,h/dTi, <Mx>= lOOGeV (b) dNj,h/dTi, <Mx> = 140 GeV 




Figure 20: Pseudorapidity density dNch/drj (in the total cms) in single diffraction in two mass 
bins with (Mx) = 100 GeV and (Mx) = 140, in pp coUisions at W = 546 GeV. The resuhs from 
DIPSY are compared with data from UA4 [3]. The excess of particles at large r] is related to the 
lack of a leading baryon in our model, see the main text. 



larger Mx-values, in fair agreement with the experiment. We also note that the simulation 
gives a surplus of particles at high rj, amounting to approximately one extra particle. This 
fact is related to the lack of a leading baryon in our model. This baryon is expected to take 
a significant fraction of the forward energy, and thus reduce the particle density for large 
r^-values. The lacking baryon is a consequence of our simple proton wavefunction, which 
contains only gluons. This was motivated by the fact that after a long cascade evolution, 
the resulting parton distribution at small x is rather insensitive to the initial wavefunction. 
The particles at high r] in fig. do, however, depend sensitively on the large x partons. 

The lack of quarks can also explain our general overestimate of the charged multiplicity, 
which is almost constant over a wide range of Mx, as can be seen in fig. 21. Our diffractive 
final states will always contain a double string connecting a leading gluon in the backward 
direction with the remnant, while in reality there will also be diffractive states with a 
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Figure 21: The average charged multiplicity as a function of Mx generated by DIPSY and com- 
pared to data from UA4 [3]. 

leading quark, giving rise to only one string and, therefore, approximately only half the 
multiplicity. The fraction of such states should decrease with increasing Mx and the rise of 
the multiplicity is mainly determined by gluonic components, but it may very well explain 
the constant surplus multiplicity in our model. In addition, the fact that a leading quark 
has a harder fragmentation function may partly explain why our simulations undershoots 
data at large negative rapidities in fig. 

The lack of quarks has also effects in other places, where large x-values are important, 
e.g. for jets with very high p±, and in future improvements we will have to include quarks 
in the proton wavefunction. 

7. Future developments 

We here discuss some effects and reactions, which are not implemented in the present MC, 
but which can be included in future versions. 

7.1 Double diffraction 

The results for single diffractive excitation, discussed in the previous sections, can also be 
generalized to double diffraction. The result is, however, numerically more complicated 
and the implementation in a MC would be more time consuming. For simplicity we here 
discuss the result for independent gluon cascades, but the result can be directly generalized 
to dipole cascades. Consider a diffractive reaction in which the projectile is excited to a real 
cascade -Rp, and the target to a real cascade i?f The number of "last" gluons is Up and rit 
respectively. The last gluons can be emitted either before or after the interaction, which im- 
plies that there are 2"'?"''"'* different diagrams contributing to the amplitude. For a specific 
diagram the sets of last gluons emitted before the interaction in the projectile and target 
cascades are denoted Bp and i?f Each pair of two last gluons in these sets, bp from the 
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projectile and bt from the target, contributes a weight Ylv,, '^Vt X^vit '^Va ^-'^P(~-^V'i,p,V'i,t)- 
Here the sums run over all virtual cascades emitted from bp and bt . (We here use the nota- 



tion introduced in sec. 4.4.3, and let Fv^^py^^ include contributions from the parent gluons 
bp and 6f ) In addition these gluons also interact with all "parent" gluons coming from the 
other side. The "last" gluons, which are emitted after the interaction, do not take part in 
the interaction, and give only the factor (—1)^, where N is the number of such gluons in 
the diagram. 

The parent gluons also interact with each other in the same way as for single diffraction. 
In obvious notation the total contribution from this diagram is 



(-1) 



N 



n ^^p"^p n ^^'"^t 

,&Rp J KneRt / 



X n E n E-^k n e*, n e 

Summing over the 2"'p different diagrams then gives the total amplitude for the double 
diffraction final state R. The result does not simplify to a factorized expression as the 
result for single diffraction in eq. ( 4.22| ), but it can still be calculated in a MC, only more 



time consuming. 

7.2 Full account of saturation effects 

In sec. 1^ the cascade evolution was assumed to be linear, that is, each gluon will radiate (or 
each dipole split) independently from the rest of the cascade. This is a good approximation 
for 7*p, with an excited photon with moderate or high , but in the case of pp scattering, 
non- linear effects play a very important role in diffractive scattering [22] . Diffractive excita- 
tion is very much suppressed in central collisions, where diffraction is dominated by elastic 
scattering. Thus diffractive excitation is largest in peripheral collisions, where saturation 
is not equally important. In the present MC, saturation effects are, besides from multiple 
subcollisions in the Lorentz frame used, also included with swings within individual sub- 
cascades, as described in sec. ^. We expect this to account for most of the saturation effects 
in events with not too high masses Mx- However, for excitation to higher masses swings 
between dipoles in different cascades may also be important, and we here discuss how it is 
possible to include these interactions, and thus take full account of saturation effects. 

This problem is similar to the problem with double diffraction discussed above, in 
that the sum of contributions from all possible diagrams does not factorize. Consider a 
real cascade with n > 2 branches ending in childless gluons U. As described in sec. 
these gluons can be emitted either before or after the interaction. This implies that there 
are 2" different diagrams which contribute to the amplitude. Virtual sub-cascades can 
be emitted from a childless gluon if this is emitted before interaction, but not if it is 
emitted after the interaction. If the sub-cascades do not interact with each other, the 
result factorizes as shown in eq. ( [4.22| ) or ( 4.24| ). This is no longer the case if the evolution 
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of one virtual sub-cascade depends on whether another childless gluon is emitted before or 
after the interaction with the target. This implies that all the 2" different diagrams must 
be calculated separately, and added to give the full production amplitude. This is not a 
problem in principle, but it becomes numerically difficult as the virtual cascades have to 
be calculated for each of the 2" possible combinations of childless gluons emitted before 
the interaction. It becomes very time consuming, unless an efficient weighting system 
is introduced, which favours the generation of states with a relatively high production 
probability. 

7.3 Quarks in the proton wavefunction 

The model for a proton wavefunction in the MC contains only gluons. The motivation is 
that gluons dominate the cascades at high energy, and for evolutions down to small x, the 
result is rather insensitive to the exact starting configuration. For diffractive events this 
motivation is only relevant for excitation to large masses. For lower Mx, experimental 
results are consistent with particle production from a single string, stretched between a 
pulled out quark and the proton remnant [30]. To describe these states it is necessary to 
include quarks in the initial proton wavefunction. This problem will be addressed in future 
work. 

7.4 Formalism combining diffractive and non-difFractive events 

We have in this paper discussed diffractive events with one or more pomerons exchanged. 
The pomerons are ladders formed by gluon pairs, which mediate momentum exchange 
and thus are able to excite the projectile to higher mass. In our earlier treatment of 
non-diffractive reactions [47], we have studied events where one or more dipole branches 
interact with the target via gluon exchange. This exchange causes a colour connection 
between the projectile cascade and the target, and provides momentum which puts the 
interacting dipole branch on shell. Branches which do not interact in this way are treated 
as virtual, and reabsorbed in the projectile state. 

We here discuss how it is possible to include exchange of both (uncut) pomerons and 
gluons in a unified formalism. An example of an event with both types of exchange in a 



single event, is shown in fig. |2^. There is no gap in the event, but the branch connected 
to the pomeron is coming on shell, increasing the multiplicity. In the formalism in ref. [47] 
the probability for this final state is included in the non-diffractive cross section, but the 
branch interacting with the pomeron is treated as virtual, and neglected in the final state. 
Although not modifying the inclusive non-diffractive cross section, this implies a small 
underestimate of the multiplicity in these events. It would, however, not affect any results 
presented in this paper. A further discussion about this combined formalism, and how it 
could be implemented in the event generator, is presented in appendix B. 

7.5 Hard diffraction and multiple gap events 

The present version of the DIPSY MC generates unbiased events. This implies that it is 
quite inefficient for producing hard diffraction, like events with high p± jets, for which some 
kind of theoretical trigger would be needed. It also does not include multiple gap events. 
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Figure 22: A cascade with botli colourful and colour neutral exchange with the target. The 
spiral line represents exchange of a single gluon with the target, which gives a colour connection. 
The zigzag line represents exchange of a colour neutral pomeron. It transfers momentum implying 
that the branch ending in the gluon marked I, can come on shell, but it does not result in a colour 
connection between this gluon and the target. 



Here we believe that the formalism described in sec. 7.4, can be generalized to describe 
events with several uncut pomerons, including difFractive events with more than one gap, 
or double diffraction with overlapping diffractive systems. 



8. Conclusions 

In this paper we describe a formalism, in which exclusive final states in difFractive excitation 
can be determined from basic QCD dynamics, BFKL evolution and saturation. Although 
diffraction is basically a quantum-mechanical phenomenon with strong interference effects, 
we show here how it is still possible to calculate the different interfering components to the 
amplitude in a MC simulation, add them with proper signs, and thus calculate the reaction 
cross section for exclusive difFractive final states. 

Our formalism is based on the Good- Walker formalism for diffractive excitation, and it 
is assumed that the virtual parton cascades represent the difFractive eigenstates defined by 
a definite absorption amplitude, in analogy with refs. [18,20-22]. Thus we here assume that 
events with large rapidity gaps can be understood as analogous to diffraction in optics, with 
matrix elements determined by absorption into inelastic channels via the optical theorem. 
This is also the case for the Regge Formalism, where the calculation oF difFractive cross 
sections is based on the AGK cutting rules and Mueller's triple-Regge formalism. In ref. 
[23] we also conjecture that the Good-Walker and triple-pomeron formalisms actually are 
different views oF the same dynamical phenomenon. 

Also a description where gaps are the result oF "soFt colour reconnection" in initially 
colour connected inelastic events, has successfully described many experimental observa- 
tions [24,25]. This is also the case For a description of gap events in terms oF projectile- 
pomeron collisions [41], which is implemented in a number oF event generators. The true 
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nature of rapidity gap events is therefore still not fully revealed. An essential feature of our 
scheme is here that the result is fully determined by the inelastic reactions, with no extra 
tunable parameters, pomeron flux factors or pomeron parton distributions, which have to 
be tuned to data. 

The formalism is based on the Lund dipole cascade model, which in turn is based on 
BFKL evolution and saturation, and is implemented in the DIPSY event generator [46,47]. 
The model is a generalization of Mueller's dipole cascade model [49-51], including also 
non-leading effects and saturation effects within the evolution. It has previously been suc- 
cessfully applied to inclusive diffractive cross sections in DIS and pp collisions [21,22], and 
later also to exclusive final states in non-diffractive events [47]. The assumption that the 
interaction is dominated by absorption into inelastic channels implies that all contribu- 
tions to the amplitudes are real, and therefore can be added with their relative signs, and 
afterwards squared to give the relevant reaction cross sections. This feature also implies 
that it is possible to calculate the Fourier transform of the 6-distribution, and determine 
the i-dependence (some results are presented in ref. [22]). 

We here note that, as discussed in refs. [60,61], the structure of non-diffractive final 
states is determined by the "backbone" or "A;j_-changing" emissions. For exclusive final 
states, final state radiation has to be added with Sudakov form factors, thus not changing 
the inclusive cross sections. Assuming that diffraction is the shadow of absorption to 
inelastic states, we here implicitly assume that also diffractive states are specified by their 
/cj^-changing backbone gluons, where final state radiation and hadronization must be added 
to get the detailed final states. 

We show in this paper some early results for single diffraction in DIS and pp collisions, 
with comparisons to experimental data. As emphasized above, these results were obtained 
without introducing any new parameters to the model, and the predictions for exclusive 
diffractive observables are entirely determined by tuning to non-diffractive and elastic data. 

In the near future we plan to present more detailed analyses of the model, include 
quarks in the proton wavefunction, and make predictions for LHC energies. We also want to 
compare our results with those from other approaches, e.g. Pythia [45] and SHRIMPS [16]. 
In the somewhat longer perspective we may include double diffraction, multiple gap events, 
and a formalism combining diffractive and non-diffractive events, as described in sec. |^. 
This development will need a significant improvement of the algorithms. Hard diffraction 
would in addition need some kind of theoretical trigger. 
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Appendices 

A. Final states for a full continuous cascade 

We here derive the results for the transition amplitudes in eqs. (|]2|,|^). For simplicity 
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Figure 23: The final state R = PUl gets contributions from diagrams where I is emitted before 
interaction (left) and after interaction (right). Due to ordering, virtual cascades from I (such as Vi) 
are only possible when I is emitted before interaction. 



we only discuss the independent gluon cascade. The modifications in a dipole cascade were 
discussed in sec. 4.4.3| . 



We first study the simpler case with a single real chain, where the set of "last" gluons, 
L, contains only a single gluon I. This is the case illustrated in fig. The emission 
structure of the real gluons R can then only be in two configurations: either they are all 



emitted before interaction (left diagram of fig. 23), or the last gluon, I, is emitted afterwards 



(right diagram of fig. ^3|). As discussed above, the "parent" gluons must always be emitted 
before interaction, to give the last gluon a chance to be emitted. 

We begin with the diagram on the left of fig. where the last real gluon is emitted 
before the interaction. This part of the 5-matrix must include a sum over all possible virtual 
cascades Vi, emitted from the partons, i, in the real cascade. The emission probability for 



a particular virtual cascade Vi is given by Cy-. defined in eq. ( 4.20 ), which satisfy the 
constraint "^Cy, = 1 where the sum runs over all virtual cascades for a fixed real parton i. 
The contribution to the ^-matrix from all diagrams of this type is given by 



5'b = /JrOr ^ 





Here fr and Fy^ = ^ fv denote the Born-level interactions, with the sum running over all 
gluons in the cascade Vr- Thus the factor e~^f'"^^^^^ represents the weight for no inelastic 
interaction by the gluon r, or by its associated virtual cascade Vr- The sums in eq. ( [A.l ) 



run over all possible virtual cascades, with their respective probabilities. In the last line 
we have separated the product over real gluons into a product over parent gluons times a 
contribution from the last gluon /. We have also pulled out the product over jdrOtr- 

The second contribution to the 5-matrix, where I is emitted after the interaction (see 



right diagram of fig. 23), contains an extra (—1), and allows only virtual cascades from the 



-33- 



"parent" gluons p. Thus this contribution is 



Sa = (-1) f n ^ra] {\{Y.-\-''''"A ■ 
\re-R / \peP Vp I 



Adding the contributions in eqs. ( A.l , |A.2| ) gives the S'-matrix element for transition 



from an incoming state |0) to an excited state \R). For a non-elastic transition we have 
T = —S, which gives the result 

n\0)^\R)) = -{Sb + Sa)= [Uf^rar] (lIE^vs, 

\reR J VpeP Vp 





n E 1 I y: cl, (1 - e-/-^-. ) I . (A.3) 

pe-P Vp 

In the last equality we have here used the relation Cy-^ = 1. 

In a general real cascade there may be several branches, which end in a "last" gluon. 
As before, the parent gluons in the real cascade, which are not the last ones in their branch, 
have to be emitted before the interaction. Otherwise their children could not be emitted. 
As in the simpler example above, they must all contribute a factor {Ylvp'^Vp^~^''~^^'')- 
The last gluons can all be emitted either before or after the interaction, and they will 
therefore all give a factor (— ^vS^ ~ ^~'^'~^'^')) to the S'-matrix element. As the sign 
of the transition amplitude is not important for the reaction rate, we get for a cascade of 
independent gluons the final result (c/ eq. ( [4.22 )) 



|T(|0) ^ \R))\ = n f^rar UE'l^'''-''^' UE-vM - ^-''-'^^) 1 (A.4) 
Vrei? / \peP Vp J \«gL Vi 

B. Events with both diffractive and non-diffractive subcoUisions 

In the Good- Walker formalism only diffractive states are considered. Cascades which get 
colour-connected to the target are treated as absorbed, and removed from the incoming 
wave. The coherent sum of cascades forming the proton wavefunction is distorted when 
different components are absorbed with different probabilities. This distortion, which leads 
to diffractive excitation, is described by a non-unitary 5-matrix, which is insensitive to the 
nature of the absorbed states. 

Within the Lund dipole cascade formalism it is also possible to describe the exclusive 
non-diffractive states responsible for the absorption [47], and it is possible to combine the 
two descriptions in a unified formalism, which can be implemented in the event generator 
DIPSY. In such a scheme the interaction will be described by a unitary S'-matrix. 

Toy model To illustrate the method we first study the simple toy model discussed in 



sec. 4.2.1. In this model there is an initial valence particle, which is able to emit a gluon. 



There are two different diffractive states, |1, 0) and |1, 1), where the numbers and 1 denote 
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an empty or occupied parton state. We now extend the set of states, and let the number 2 
indicate a parton which has interacted via gluon exchange, and is attached to the target. 
Thus we have 6 different states: 



|1,0), |1,1), |1,2), |2,0), |2,1), |2,2). 



(B.l) 



We have also to generalize the evolution operator U{—6,—oo) in eq. (4.3), which de- 
scribes the emission of a gluon. The valence particle can emit a gluon with weight /3, if the 
gluon state is empty. This is possible also if the valence parton is attached to the target, 
and we therefore have 



C/(-<5,-oo)|l,0) =q|1,0) +/3|1,1), f/(-5,-oo)|2,0) = a|2, 0) + /3|2, 1). 



(B.2) 



All other states already contain an emitted gluon, and are unaffected by U{—5, -co). 

The interaction with the target described by U\at has now, besides the elastic non- 
interaction weights e~-^° and e~^^ in eq. (4.4), also non-diffractive interactions transforming 



a free valence parton to an attached one, with weight ^/l — 6'"^^" , and similarly transform- 
ing a free gluon to an attached gluon with weight Vl — e~2/i . This implies that Ui^t can 
be written in a factorized form U{—5, +5) = Umt = fAnt,of^int,i, where 



f^int.oll,^) =7o|l,^) +'5o|2,i), f^int,l|j, 1) =7i|i, 1) +'5i|j,2) 



(B.3) 



Here 7i = e -^^ is the non-absorption amplitude, and 6i = \/l — e "^f' the non-diffractive 
interaction amplitude. C/int,o does not act on gluon 1 and vice versa. Note that 

l! + 6f = l, 

as now C/int is a unitary operator. 

Finally the evolution after the interaction is given by U{+oo, +6) = W {—6, —oo), 
which implies that the non-diagonal elements have changed sign. (We have defined the 
phases such that a and /3 are real.) Note that since U{—6, —oo) can only emit free gluons, 
an attached gluon cannot be absorbed by U{+oo,+5). The unitary S'-matrix is given 
by S" = {—6, —oo)UintU{—S, —oo), and it is straightforward to read off the transition 
amplitudes for an incoming state with a single valence parton: 



S(|1,0) - 


Hi,o)) 


S(|1,0)- 




S(|1,0)- 


>|1,2)) 


5(|1,0)- 


> |2,0)) 


S(|1,0)- 


>|2,1)) 


5(|1,0)- 


> |2,2)) 



7o(a2 _^ ^2^_^^) ^ g-/o ^^2 _^ ^2g-/i^ 

— q/37o(1 — 7i) = —af3e~^° ^ 



1 - e 



-a/35o{l - 71) = -a/3\/l - e-2/0 (l - e'f^ 
l35o5i = /3Vl- e-2/0 Vl-e-2/1 
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We note here that the amphtudes for elastic scattering and diffractive excitation, in the 
first two lines, agree with the result in sec. |4.2.1 . The inclusive non-diffractive cross section 



also agrees with the eikonal result dai^ei/d b = (1— e" 



-2F\ 



-2/0 



-2/0-2/1 



used in earlier publications. However, in our analysis of exclusive non-diffractive final states 
in ref . [47] , the possibility for momentum exchange which can put a virtual branch on shell, 
as the emitted gluon in state |2, 1), was not included. Thus in the present version of the 
DIPSY event generator the cross section for the state |2,0), with only the valence parton 
connected to the target, is overestimated and given a probability representing both states 
|2, 1) and |2,0). 

General cascade The above result generalises easily to a full cascade of multiple chains 
with continuous emissions. For clarity we here show the result for cascades of independent 
gluons, as discussed in sec. |4.3| . We let denote the non-diffractive cascade, consisting of 
the interacting gluons (set /) and their ancestors (set A). As before R is the set of gluons 
in the branches which come on shell via pomeron exchange, with L denoting the last gluons 
in its branch and P denoting their parents (including grandparents etc., but not including 
any gluon in N) . The generalization of eq. ( |4.22[ ) then reads 



5(|0) ^\N + R)) 




\aeA Vo 



vrei? / \p6P Vp J \ieL V, J ) 

This result can be implemented in a future version of DIPSY, and would imply slightly 
higher multiplicities in non-diffractive events. 
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